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Abstract 

The quantum kicked particle in a magnetic field is studied in a weak-chaos 
regime under realistic conditions, i.e., for general values of the conserved coor- 
dinate x c of the cyclotron orbit center. The system exhibits spectral structures 
["Hofstadter butterflies" (HBs)] and quantum diffusion depending sensitively 
on x c . Most significant changes take place when x c approaches the value 
at which quantum antiresonance (exactly periodic recurrences) can occur: 
the HB essentially "doubles" and the quantum-diffusion coefficient D{x c ) is 
strongly reduced. An explanation of these phenomena, including an approx- 
imate formula for D(x c ) in a class of wave packets, is given on the basis of 
an effective Hamiltonian which is derived as a power expansion in a small 
parameter. The global quantum diffusion of a two-dimensional wave packet 
for all x c is briefly considered. 
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I. INTRODUCTION 



Simple quantum systems whose classical counterparts are nonintegrable exhibit a variety 
of remarkable phenomena [1] which continue to attract much interest in several contexts. A 
realistic system, whose experimental realization by different methods has been the object of 
many recent works (see a very partial list in Ref. [2]), is the well-known kicked rotor [1,3,4], 
modelling the dynamical-localization phenomenon. A class of systems exhibiting basically 
different phenomena, such as quantum diffusion associated with a fractal spectrum, are rep- 
resented by Hamiltonians periodic in phase space [5-20]. The simplest such system is the 
integrable Harper model [5] with Hamiltonian H = cos(:r) + cos(p), whose quantization de- 
scribes the energy spectrum within a Bloch band in a magnetic field or within a broadened 
Landau level in a crystalline periodic potential [5-7,10]. Plotting this spectrum as a function 
of the magnetic flux through a unit cell (this flux is related to a scaled Planck constant for 
the Harper model) gives the well-known "Hofstadter butterfly" [6], featuring a clearly fractal 
structure [9]. Quite recently [10], much of this structure was detected experimentally. 

The simplest nonintegrable system with a mixed phase space and a periodic Hamiltonian 
is the kicked Harper model (KHM) [11-19] with H = cos(p)+cos(x) J2T=-oo ^/ T - s )> whicn 
reduces to the Harper Hamiltonian when the time period T — > 0. The KHM models the 
regular-chaotic transition [11,19] as well as the effect of chaos on a fractal spectrum and on 
quantum diffusion [13]. It seems that the only realistic interpretation of the KHM is its exact 
relation [16] with a particular case of the two-dimensional (2D) system of the periodically 
kicked particle in a magnetic field [21-23]. This system is described by 



where II = (U x , U y ) = p— B x r/ (2c) is the kinetic momentum of a particle with unit charge 
and unit mass in a uniform magnetic field B (along the z-axis) and K is a nonintegrability 
parameter. Let us first summarize some known facts about (1) [21-23]. It is crucial to 
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represent (1) in the natural degrees of freedom in a magnetic field, given by the independent 
conjugate pairs (U x , U y ) and (x c , y c ) (coordinates of the center of a cyclotron orbit) [24]. 
Defining u = H x /u, v = H y /u, where uj = B/c is the cyclotron frequency, and using the 
relation x c = x + H v /uj = x + v (from simple geometry), (1) can be rewritten as follows [22] 



Since y c does not appear in (2), its conjugate mate x c is conserved and can therefore be 
treated as a parameter. This reduces (1) to an ensemble of periodically kicked harmonic 
oscillators parametrized by x c . Classically, these degenerate systems are expected to exhibit 
unbounded chaotic diffusion in the (u, v) phase plane for arbitrarily small values of K, 
especially under resonance conditions, with a = uT a rational multiple of 2n [21,22]; this 
diffusion is observed to take place on a "stochastic web" whose symmetry depends on a. Ex- 
amples of square-symmetry webs (a = rr/2) are shown Fig. 1. The KHM is exactly related 
to (2), in essence, only for a = n/2 and for very special values of x c , x c = 0, n (see Sec. A of 
Appendix and generalizations of this relation in Refs. [15,16]). Now, the classical properties 
of (2), in particular the structure of the stochastic web, its width, and the diffusion rate on 
it, are known to depend strongly on x c [22,23] (see also Sec. III). In addition, one should 
notice that realistic classical ensembles or 2D wave packets for (1) exhibit all values of x c . 
These facts motivate the investigation of the quantum properties of (2) for arbitrary x c . 
Except of general results [15,16], these properties appear to be essentially unexplored. 

In this paper, we present a first study of the quantized system (2) for general x c . We 
focus on the weak-chaos regime of small K and, as in the classical studies [22,23], on the 
case of a = ir/2. The global spectral features of the system at fixed x c are exhibited by a 
suitably defined plot of the quasienergy spectrum as a function of a scaled Planck constant 
h. The term "Hofstadter butterfly" (HB) is extended so as to refer to such a plot. The quan- 
tum properties are found to depend sensitively on x c . Most significant changes take place 
when x c approaches the value of 7i/2 at which quantum antiresonance (QAR), i.e., exactly 
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periodic recurrences [15,25], occurs for integer h/(2n). While the HB for most values of x c 
is approximately the standard one of the Harper model [6], the HB for x c = it/ 2 is, up to 
a local scaling, a perturbed "doubled" version of the standard HB. This change of the HB 
structure in the "QAR transition" x c — > it/ 2 is explained on the basis of an effective Hamil- 
tonian which is derived as a power expansion in a small parameter. For irrational h/(2n), 
the time evolution of wave packets is observed to exhibit approximately an asymptotic quan- 
tum diffusion. This is characterized by either a diffusion coefficient D(x c ) at fixed x c or a 
global coefficient D, given by a weighted average of D(x c ) over x c . An approximate formula 
for D(x c ) in a simple class of wave packets is derived from the effective Hamiltonian and is 
verified numerically. As x c — > tt/2, D(x c ) decreases and QAR is manifested by the fact that 
D{tt/2) is significantly smaller than typical values of D(x c ). While these phenomena are of 
a purely quantum nature, they have classical analogs. The paper is organized as follows. In 
Sec. II, we recall the QAR phenomenon [15,25] and study general spectral properties of the 
system (2), also on the basis of an effective Hamiltonian which is derived as an expansion in 
a small parameter. In Sec. Ill, we consider classical limits and analogs of some of the results 
in Sec. II. In Sec. IV, we study the quantum diffusion of wave packets for general x c , using 
also the effective-Hamiltonian approximation. A summary and conclusions are presented in 
Sec. V. Proofs of several statements in Sees. I, II, and IV are given in the Appendix. 



II. QUANTUM ANTIRESONANCE, GENERAL SPECTRAL PROPERTIES, 

AND EFFECTIVE HAMILTONIAN 

To quantize the system (2), (u, v) are replaced by the corresponding conjugate operators 
(w, -0). From the definitions u = H x /w, v = H y /uj, where (II X , n^) = p — B x r/(2c) and 
uj = B/c, it is easy to see that (it, v) satisfy [u, v] = ih/uj. We choose units such that uj — 1, 
so that [u, v] = ih, where K is a scaled Planck constant. The one-period evolution operator 
for the quantized system (2) will be denoted by Ut(x c ). Let us now recall the phenomenon of 
quantum antiresonance (QAR) [15]. Under classical resonance conditions, a = uT = 2irj/l 
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(j and / are coprime integers), QAR for the system (2) is defined by U l T (x c ) = constant phase 
factor, i.e., exactly periodic recurrences with the natural resonance period IT. As shown 
in Ref. [15], QAR occurs in the nonintegrable (I > 2) case of a general kicked harmonic 
oscillator only if three conditions are satisfied: (a) The kicking periodic potential is an odd 
function, up to an additive constant; this implies that x c = ir/2, 3tc/2 for the potential 
cos(x c — v) in Eq. (2). (b) I — 4 or I — 6, corresponding to stochastic webs with square 
or hexagonal symmetry, respectively, (c) h/(2n) is integer for I — 4 while \f3h/(Aic) is 
integer for I — 6. When all these conditions hold, it follows from the results in Ref. [15] 
that U l T {x c ) = —1 for x c = tt/2, 3tc/2 and I — 4, 6. In this paper, only the case of / = 4 
with j/l = 1/4 (a = 7r/2) will be considered. Using general results [15], one finds that the 
evolution operator U^(x c ) in this case is equal to —U(x c ), where 

U (x c ) = exp[i/x cos(x c — u)] exp[i/x cos(x c + -0)] exp[i/x cos(x c + u)] exp[i/x cos(x c — v)] (3) 

and fx = K/h. We show in the Appendix (Sec. A) that the range of x c in Eq. (3) can be 
restricted to [0, tc/2] without loss of generality. Thus, only the QAR point x c = n/2 will be 
considered; U(n/2) = 1 for integer h/(2n). 

We now extend the term "Hofstadter butterfly" (HB) [6] to the operator (3). Let us 
denote by E and respectively, the quasienergies (QEs) and QE eigenstates, defined by 
U(x c )^ = exp(-iLiE)^, with E lying in the interval [—tt/li, tt/li). We then define the 
HB as the plot of the QE spectrum E(h) as a function of h, for < H < 2n, at fixed \i. 
Such a plot was introduced in Ref. [13] for the KHM and is motivated by the fact, shown 
in the Appendix (Sec. B), that E(h) is 27r-periodic in h at fixed ii. A small value of il 
characterizes a HB associated with a "semiclassical weak-chaos regime", i.e., a regime of 
small K — h/i for all small values of h, h < 1. We also show in the Appendix (Sec. B) that: 
(a) For rational h/(2n) = q/p (q and p are coprime integers), the QE spectrum consists of 
p bands, each g-fold degenerate, (b) The HB has reflection symmetry around H — n [i.e., 
E(2n — h) = E(h)] only if x c = or x c = tc/2. Examples of HBs, plotted for all h = 2nq/p 

5 



with 1 < q < p < 50, are shown in Fig. 2 for three values of x c . We see that for x c = tc/2 
(Fig. 2(c)) the QE spectrum E — > as h — > 27r or h — > 0, reflecting the QAR phenomenon, 
t>(vr/2) = 1. 

The unitary operator (3) can be formally written as U(x c ) = exp[— ifj,H e s(x c )], where 
H e s(x c ) is an effective Hamiltonian for the problem. Clearly, if E' is the "energy" spectrum 
of H e g(x c ), the QE spectrum E is just E' modulo the interval [— 7r//i, 7r//i). Thus, the 
"unfolded" HB is the spectrum E' plotted as a function of h (0 < h < 2ir) at fixed //. The 
dependence of some general HB features on x c can be understood on the basis of an expres- 
sion for H cS (x c ) which we now derive. We first note that the arguments of the exponents in 
Eq. (3) involve simple operators exp(±ra) and exp(±it>), giving just translations by ±h in 
the (u, v) phase plane. Using then the known formula [26] 

exp(A) exp(5) = exp (a + B + ±[A, B] + ±[A, [A, B]] + ±[[A, B], B] + ..?j, (4) 

valid for arbitrary operators A and B, one can verify by tedious but straightforward algebra 
that H cS (x c ) can be expressed as a power expansion in a parameter e, 

H cS (x c ) = fytf r (:r c ), 6 = vi8m{h/2) = f S ^01. ( 5 ) 

r=0 1 

The zero-order coefficient in Eq. (5) is 

Hq(x c ) = — 2 cos(x c )[cos('u) + cos(t))]. (6) 

The operator (6) is a Harper Hamiltonian and if e is sufficiently small, 2| cos(x c )/e| ^> f, 
is the dominant term in the expansion (5) for most values of x c . For general /x, e is small 
provided h/ (2tt) is sufficiently close to an integer. In a weak-chaos (small fi or K) regime, 
on which we shall focus, e is globally small for all H. Then, for 2| cos(x c )/e| 3> 1, the HB is 
essentially the standard one for the Harper model [6], see Fig. 2(a). For 2| cos(a; c )/e| < f, 
higher-order coefficients in Eq. (5) become significant, especially for h ~ tt (largest e at 
fixed /i), leading to a clearly non-standard HB with a visible reflection-symmetry breaking; 
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see Fig. 2(b). 



For x c = 7r/2, H vanishes. If also h/(2ir) is integer, e = and thus H e s = 0, implying 
U(n/2) = 1 (QAR). Small values of e for x c = n/2 define a QAR vicinity; thus, a weak- 
chaos regime at x c = n/2 corresponds to a QAR vicinity for arbitrary h. Now, H cS (n/2) = 
eH\ + e 2 H 2 + • • •, and we find that 

Hi(n/2) = -[cos(u + v) + cos(w - v)]. (7) 

Since the operators v! = u + v and v' — v — u are clearly conjugate, [u', v'} = ih' with 
h' = 2h, we see that Hi in Eq. (7) is again a Harper Hamiltonian. When h varies from to 
27r, h! = 2h varies from to Atx. Then, the HB for H eS /e — Hi + eH 2 + ... is a perturbed 
doubled version of the standard HB, see Fig. 3. A simple measure of the perturbation of 
this HB is P = \l-AE'(n)/AE'(2n)\, where AE'(h) is the width of the spectrum of H eS /e. 
The exact result P = 1 — fi~ 2 arcsin[sin 2 (/i)] = /i 2 /3 + . . . is easily derived (see note [27]). As 
h approaches its QAR value of 2ir (or at fixed /i), the width AE(h) ~ /ism(h/2)AE'(2ir) 
of the QE spectrum vanishes almost linearly in h, see Fig. 2(c). 

A simple consequence of the approximately doubled structure of the HB for H e ^{n/2)/e 
is that if h = 2nq/p and p is even, p = 2p', the p bands in this HB actually form p' pairs 
of overlapping bands. This is because the spectrum of the Harper Hamiltonian (7) for 
h' = 2h = 2-Kq/p' always consists of p' nonoverlapping bands [6]. Each such band must then 
correspond to a pair of overlapping bands of H cS for sufficiently small e. Thus, the number 
of gaps in the spectrum for even p is more than halved when x c is varied from to ir/2. 

III. CLASSICAL LIMITS AND ANALOGS 

In this section, we consider classical limits and analogs of some of the results above. The 
classical one-period map for the system (2) in the case of a = n/2 is [22]: 
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M : u s+ i = v s , v s+ i = -u s + K sin(:r c - v s ), (8) 

The fourth iterate M 4 of M is the classical map corresponding to the evolution operator 
U(x c ) in Eq. (3). One can easily show that M 4 is a near-identity map for small K [23], i.e., 
u t+ i — u t + O(X), f t+ i = v t + O(i^), where t denotes now an integer time index counting 
iterations of M 4 . It is then clear from the expression U(x c ) = exp[—iKH cS (x c )/h] that the 
classical limit (h — > 0, e — > if/2) of H cS (x c ) is an integrable Hamiltonian H e g(u, v; x c ) 
generating, in a continuous time denoted here by t', a phase-space flow which approximates 
the near-integrable map M 4 at times t' = Kt: 

dH cS _ du ^ u t+1 - u t dH eS _ dv_ ^ v t +i - Vt / g x 

Expressions for KH e g/A as power expansions in K are given in Ref. [23]. Like U(x c ), M 4 
and H e ff(u, v; x c ) are generally periodic in the phase space (u, v) with a unit cell of area 
An 2 , so that the orbit structure of M 4 or the flow (9) feature the same periodicity, as clearly 
shown in Figs. 1(a) and 1(b). For x c = n/2, however, the dominant term in H cS (u, v; vr/2) 
is the classical analog Hi of (7) which is periodic with a unit cell of area 2tt 2 , half the size 
of the usual unit cell; this periodicity is then approximately exhibited by the corresponding 
orbit structure, compare Fig. 1(c) with Fig. 1(a). This is the classical fingerprint in the 
approximate doubling of the HB for x c = n/2 (see previous section and Fig. 3). 

The hyperbolic fixed points of M 4 form a periodic array connected by the stochastic 
web, see Fig. 1. As shown in Ref. [23], the width of the web for small K is proportional to 
exp[— 7r 2 / ln(A)], where A is the largest eigenvalue of the linearization of M 4 at a hyperbolic 
fixed point and is given by A 1 + 2ifcos(a; c ) for cos(a; c ) 3> K/A and by A 1 + K 2 for 
x c = 7r/2. Thus, the latter case is characterized by a significantly small web width and slow 
chaotic diffusion relative to those for most values of x c [23]. These features may be viewed 
as classical analogs of the QAR phenomenon, including some quantum behaviors in a QAR 
vicinity (see Sec. IV C and Sec. V). 
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IV. QUANTUM DIFFUSION FOR GENERAL x c 

In this section, we study the quantum diffusion exhibited by wave packets evolving under 
the operator (3) and its effective-Hamiltonian approximations. 

A. General quantum evolution and periodic wave packets 

A general wave packet |<3>) for the original 2D problem (1) can be naturally expressed in 
a representation based on the two degrees of freedom (conjugate pairs) (u, v) and (x c , y c ), 
for example, in the (u, x c ) representation, (u, x c \$) = x c ). The time evolution of 

$(«, x c ) can be decomposed into the independent evolutions for the separate, conserved x c 
values, 

$(u, x c - t) = U\x c )^(u, x c ), (10) 

t being the integer time index. We show in the Appendix (Sec. C) that the quantum 
dynamics (10) can be fully reproduced from that of wave packets on a cylindrical phase 
space, < u < 2n, — oo < v < oo, with <£>(«, x c ) 27r-periodic in u: 

oo 

$(«, x c ) — E] (p(n, x c )exp(inu). (11) 

n=— oo 

Here the Fourier coefficients ip(n, x c ) give the (v, x c ) representation of |$) on the cylinder, 
with v quantized in units of H, v — nh. The quantum evolution of the wave packet (11) in 
the (n, x c ) representation, 

ip(n, x c ; t) = C/*(x c )</?(n, x c ), (12) 

can be easily calculated as in the kicked- rotor case [3], by using a fast Fourier transform 
to switch from the position (u) to the momentum (v = nh) representation (and vice-versa) 
before the application of each factor in Eq. (3). The spreading of the wave packet (12) along 
the cylinder at fixed x c is measured by the expectation value (v 2 ) t of v 2 in ip(n, x c ; t): 

-2 oo oo 

(v 2 ) t = E n "\^ x ^ ^l 2 ' ^c) = E ^c)| 2 . (13) 



The expectation value of v 2 in the original 2D wave packet ip(n, x c ; t) (both n and x c are 
variables) is given by 

oo „ 

1^y t = tf J2 n 2 \<p(n, x c ; t)\ 2 dx c , (14) 

n=— oo 

where normalization of tp over (n, x c ) is assumed, J N(x c )dx c = 1. We also show in 
the Appendix (Sec. C) that our main results below concerning (v 2 ) t can be extended to 
the expectation value {u 2 (u 2 + v 2 ) / 2) t of the kinetic energy in an arbitrary wave packet 
(10). Thus, for convenience and without loss of generality, we shall restrict our attention to 
periodic wave packets (11) and to the expectation values (13) and (14). 

B. Quantum diffusion at fixed x c and global quantum diffusion 

We have seen in Sec. II that if e is sufficiently small the dominant term in the expan- 
sion (5) of H e s(x c ) for most values of x c , including x c = n/2, is a Harper Hamiltonian. It is 
known [8,9] that for irrational hj {2n) the latter Hamiltonian exhibits a fractal spectrum and, 
asymptotically in time, an approximate quantum diffusion of the second moment of a wave 
packet. We can therefore expect a similar behavior of (v 2 ) t under U(x c ) = exp[— ifiH c s(x c )} 
for irrational h/(2it) and sufficiently large t: (v 2 ) t 2D(x c )t, where D(x c ) is the diffusion 
coefficient at fixed x c . 

As an illustration, let us choose the initial wave packet (11) as a periodized coherent 
state centered on some x c -dependent point (u, v), 

, . ,,1/4 ^ f ,v(x c )(u — 2irm) \u — 2nm — u(x c )] 2 1 , . 

x c ) = (nh) 1/4 ex P { 1 fi " " 2 ft } • ( 15 ) 

JTt= — OO ^ 

For sufficiently small h, a classical quantity analogous to the expectation value (v 2 ) t in the 
wave packet (15) is the average {v 2 ) of v 2 over an ensemble of initial conditions (wo? fo) 
uniformly distributed in a disk centered on (u, v) and of radius V2h; (u t , v t ) is the tth 
iterate of (uq, vq) under the map M 4 (see previous section). Fig. 4 shows log- log plots of 
(v 2 ) t and (v 2 ) for h/(2n) = [51 + (V5 - l)/2]~\ K = 0.157, and < t < 10 5 in the three 
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cases of x c considered in Figs. 1 and 2; in each case, (u, v) was chosen as a hyperbolic 
fixed point of the corresponding map M 4 . We see that {v 2 ) t and (v 2 ) essentially coincide 
up to some time t = t*(x c ), after which (v 2 ) stops spreading and there occurs a crossover of 
(v 2 ) t to an approximate quantum diffusion, (v 2 ) t ~ 2D(x c )t. In the case of x c = ir/2, this 
diffusion is shown more clearly by the normal plot of (v 2 ) t in the inset of Fig. 4. 

These results can be understood as follows. The irrational value of h/ (2ir) chosen above 
represents a semiclassical regime in which h = 2irh is small relative to the area 4n 2 of a unit 
cell but large relative to the width ~ exp[— n 2 / ln(A)] of the stochastic web for K = 0.157 (see 
classical details in the previous section). Thus, classical chaos leaves almost no fingerprints 
in the quantum properties which should then resemble those of the integrable approximation 
of the system given by the effective Hamiltonian H ef[ (x c ), with the stochastic web replaced 
by its separatrix "skeleton" [21,23]. A classical ensemble of width ~ V2h centered on a 
hyperbolic fixed point will lie mostly in the stable regions within the unit cells and will 
spread to its maximal extent (of the order of the size of one unit cell) in a time t = t*(x c ) 
roughly proportional to ln(l/ft)/ln[A(a; c )]; this agrees with numerical estimates of t*(x c ) in 
all cases considered, including those shown in Fig. 4. The quantum wave packet, however, 
will continue to spread due to tunnelling between neighboring unit cells, leading to quantum 
diffusion [8,9]. 

Assuming the asymptotic quantum-diffusion behavior of (v 2 ) t , (v 2 ) t ~ 2D(x c )t, to hold 
approximately for all x c , one finds from Eqs. (13) and (14) that the 2D wave packet (12) 
can be characterized by a more realistic quantity, its global diffusion coefficient D, given by 
a weighted average of D(x c ) over x c : 

l^Y t ^2Dt, D = J N(x c )D(x c )dx c , (16) 

where N(x c ) is defined in Eq. (13). The global diffusion (16) will be briefly considered at 
the end of this section. 
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C. x c -dependence of the quantum-diffusion coefficient for simple wave packets 

We now consider the ^-dependence of the quantum-diffusion coefficient D. It is clear 
from Eq. (12) that this dependence is generally due to that of both U(x c ) and the initial 
wave packet (p(n, x c ). For simplicity, we shall restrict ourselves from now on to the class of 
wave packets for which one has a separation of the (n, x c ) variables, ip(n, x c ) = x(x c )ijj(n) 
\x( x c) an d ipiji) are assumed to be normalized]. Then, <p(n, x c ; t) = x^^^^x^in), so 
that (v 2 ) t in Eq. (13) is independent of x( x c) an d the x c -dependence of (v 2 ) t for fixed ij)(n) is 
completely due to that of U (x c ). This allows one to derive an approximate formula for D(x c ) 
in the case of 2| cos(:r c )/e| 3> 1. In this case, which covers most values of x C) the effective 
Hamiltonian H c q Hq, where Ho is given by Eq. (6). The tth power of the approximate 
evolution operator U = exp(—ifiH ) can be written as follows: 

Uq(x c ) = exp {2iK cos(x c )t [cos(u) + cos(-O)] /h} . (17) 

Rel. (17) implies that the wave packet ip(n, x c ; i) = Uq(x c )iP(ti) depends on K, x c , and t only 
through the combination K cos(x c )t. Thus, if one assumes an asymptotic quantum- diffusive 
behavior, (v 2 ) t ~ 2D(x c )t, accurately described by U , D(x c ) must be approximately linear 
in K cos(x c ): 

D(x c ) « 2D ( n ) Kcos{x c ), 2| cos(:r c )/e| > 1. (18) 

Here the proportionality constants and are associated with the x c -intervals 

where cos(x c ) > and cos(x c ) < 0, respectively; assuming that K > for definite- 
ness, one has > and < since D(x c ) > 0. The constants and \D^\ 

are just the diffusion coefficients for (v 2 ) t under the standard Harper evolution opera- 
tor {7 H = exp {i [cos(w) +cos(v)]/H} and its inverse U^ 1 , respectively For general ip(n), 
7^ \D^ \ . However, if ip(n) satisfies some symmetry properties, e.g., its u representa- 
tion is given by the right-hand side of Eq. (15) with u a multiple of n (as in the example 
below), one can show that = ^| (see Appendix, Sec. D). Formula (18) can then 
be rewritten as D(x c ) ~ 2DuK\ cos(x c )| with the same constant for all x c satisfying 
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2\ cos(x c )/e| 3> 1. 

As a check of formula (18), we have calculated D(x c ) for h/(2n) = [51 + (V5-l)/2]-\ 
K = 0.157, and x c = kn/AO, k = 0, 1, . . . , 18, by a linear fit to (v 2 ) t (0 < t < 10 4 ), evolving 
under the exact operator U(x c ). The u representation of the initial wave packet ip{n) is 
given by the right-hand side of Eq. (15) with (u, v) = (ir, 0), a hyperbolic fixed point for 
x c = (see Fig. 1(a)). The results, shown in Fig. 5, agree satisfactorily with a least-square 
fit of the function (18) to the data. 

Consider now the case of x c = ir/2. From the fact that H cS (ir/2) eH 1 (ir/2) = O(K), 
one expects, using arguments similar to those leading to formula (18), that D(n/2) should be 
approximately quadratic in K. In Table I, we present results concerning the if- dependence 
of D(0)/K and D(tt/2)/K 2 . The values of D(0) and D(n/2) were calculated by a linear 
fit to (v 2 ) t (0 < t < 10 4 for x c = and < t < 10 5 for x c = tt/2), evolving under the 
exact operator (3); for both x c = and x c = n/2, the initial wave packet is given by the 
right-hand side of Eq. (15) with (w, v) chosen as a hyperbolic fixed point of the correspond- 
ing map M 4 . The results for D(0)/K agree very well with the approximate linearity in K 
predicted by formula (18). On the other hand, D(n/2) appears to decrease faster than K 2 
as K decreases, which may indicate that high-order corrections to H c q for x c = tt/2 have a 
stronger impact than for x c = 0. In any case, Table I provides evidence that as K decreases 
D(n/2) becomes significantly smaller than typical values of D(x c ). 

We end this section by a brief consideration of the global quantum diffusion (16) in 
the case of the initial wave packets y?(n, x c ) = x{ x c)ip{n). For such a wave packet, for- 
mula (18) can be used to calculate approximately the global diffusion coefficient D, with 
N(x c ) = |x(^c)| 2 in Eqs. (13) and (16). Choosing N(x c ) uniform, N(x c ) = 2/tt, for 
< x c < 7r/2 and vanishing outside this interval, one finds from Eqs. (16) and (18) that 
D « AD^K/tt. Much smaller values of D are obtained if N(x c ) is strongly localized around 
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x c = 7r/2. Fig. 6 shows the time evolution of {v 2 ) v given by a uniform average of {v 2 ) t over 
x c = kn/AO, k = 0, 1, . . . , 18, where {v 2 ) t for each of these x c values was calculated as 
specified above for obtaining the results in Fig. 5. We also show in Fig. 6 the functions 
{v 2 ) t for x c = (k = 0) and x c = 2n/5 (k = 16). As expected, averaging over x c removes 
from (v 2 ) t much of the fluctuations exhibited by (v 2 ) t at fixed x c . 

V. SUMMARY AND CONCLUSIONS 

Our first study of the quantized system (2) for general x c has revealed the significant 
impact of the transition x c — > n/2 to a QAR vicinity on the spectral and quantum-diffusion 
properties in a weak-chaos (small K) regime. In this transition, the dominant term in the 
effective-Hamiltonian expansion (5) changes from H to eH 1 in a x c -window of width ~ e 
around x c = n/2, where H vanishes. Both H (x c ) and H\{ir/2) are Harper Hamiltonians, 
but H 1 (n/2) is periodic in phase space with a unit cell half the size of that of H (x c ). As a 
consequence, the HB for x c = n/2 is essentially a perturbed doubled version of the standard 
HB for the Harper model (see Fig. 3). Thus, clS CCq, is varied from to tt/2, there must occur 
strong changes in the HB structure and the spectral properties (see Fig. 2 and end of Sec. 
II), reflecting corresponding changes in the classical web structure due to bifurcations [23] 
(see Fig. 1). A more detailed investigation of the fingerprints of these bifurcations in the 
quantum properties is planned for a future work. 

The quantum diffusion of wave packets for the system (2) was characterized by a coeffi- 
cient D(x c ) at fixed x c and by a global coefficient D. Formula (18) for D(x c ) was derived 
on the basis of the approximate effective Hamiltonian H Q (x c ) and it agrees reasonably well 
with numerical results for sufficiently small K (see Fig. 5 and Table I). This formula and 
numerical evidence indicate that D(x c ) is strongly reduced as x c — > n/2: D(n /2) appears 
to be not larger than 0(e) relative to typical values of D(x c ). This phenomenon is of a 
purely quantum nature but it is analogous to the relatively slow classical chaotic diffusion 
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for x c = 7r/2 in a weak-chaos regime [23]. An interesting but apparently difficult problem 
is to obtain a refined formula for D(x c ), extending the zero-order result (18) to the small 
:r c -window around x c = n/2 where higher-order terms in the expansion (5) must be taken 
into account. 

The strong variation of D(x c ) with x c is exhibited in practice by the global quantum 
diffusion of a 2D wave packet for all x c . If the potential cos(x) in Eq. (1) is replaced by 
cos(x — 9), where 6 is an adjustable phase, one can "filter" the wave-packet component 
associated with an arbitrary value of x c ~ 6 + n/2 as the component having the smallest 
quantum-diffusion rate. The new quantum phenomena predicted by the realistic, general-x c 
approach to (1) may be observed in possible experimental realizations of this system. 
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APPENDIX 

A. Exact relation between system (2) and KHM and relevant x c range 

The evolution operator (3) for the system (2) in the case of a = rc/2 is equivalent to 
U'(x c ) = exp [— i/iW(v)] U(x c ) exp [i//W(t))], where W(v) = — cos(x c — v). One has 

U (x c ) = f^KHM^KHM' (^9) 

where 

C/khm = ex P [—inW(±v)\ exp [—ijjiW(±uj\ 



15 



is the one-period evolution operator for the generalized KHM [16,17] with Hamiltonian 
H = W(±v) + W(±u) - s). This reduces to the usual KHM Hamiltonian, 

with W(y) = cos(f) (up to a sign), only for x c = 0, it. Then, f/^HM = ^khm = ^khm and 
the exact relation U'(x c ) = U^hm follows from Eq. (19). 

The evolution operator (3) satisfies also the similarity relations 

U(x c + tt) = DU(x c )D-\ U(-x c ) = SU(x c )S-\ (20) 

where 

D = exp(iriu/H) exp(mv/h), S = exp[— i/i cos(x c + -0)] exp[— i/i cos(:r c — u)]. 

Since U(x c ) is 27r-periodic in x c , Rels. (20) imply that the range of x c can be restricted to 
[0, 7r/2] without loss of generality. 

B. Quasienergy band spectrum and some of its basic properties 

Using Rel. (19), the results in Ref. [17] concerning the quasienergy (QE) band spectra 
of generalized KHMs can be straightforwardly extended to the operator U'(x c ), which is 
equivalent to the evolution operator U (x c ) in Eq. (3) and thus has the same QE spectra 
as those of U(x c ). For h = 2irq/p (q and p are coprime integers), U'(x c ) and the phase- 
space translations D\ = exp(2niu/h) and D 2 = exp(ipv) form a complete set of commuting 
operators. Their simultaneous eigenstates are given, in the v representation, by 

p— 1 oo 

^b,™,x c ( v ) = X^ b ( m ' w; Xc "> ex v[in(w 1 +mh)/q]5(v -w 2 + 2im/p). 

m=0 n=—oo 

Here the index b = 1,..., p labels p QE bands, w = (w±, w 2 ) is a Bloch wave vector 
spanning a band and ranging in the "Brillouin zone" < W\ < h, < w 2 < 27r/p, and 
{0b(m; w; x c )}^} , b = 1, p, are p independent vectors of coefficients. The correspond- 
ing eigenvalues of U'(x c ) are exp[— ifj,Eb(w; x c )], where Eb(w; x c ) is a QE band. It is easy 
to see that U'(x c ) commutes with D = exp(27ri-0 / H) and also that Dq = D 2 . Then, the q 
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states ^^j. = DQ^b,w,x c , d — 0, g — 1, are different eigenstates of U'(x c ), all associated 
with the same eigenvalue. Each band is therefore g-fold degenerate. 

The p vectors {exp(imw 2 )0fe(m; w; x c )}^ , 6 = 1,..., p, are the eigenvectors of a p x p 
w-dependent matrix M(w; x c ), 

M(w; x c ) = M (+) (w; x c )M ( -\w; x c ), (21) 

where M^fw; x c ) are matrices corresponding to the operators C/^hm i n R e l- (19) an d 
whose elements can be explicitly written using the results in Ref. [17]: 

oo 

M 2L( w ; x o) = exp[ificos{x c =f w x =fm'h)] J^ gp +m'-m)(n, x c ) exp(igpw 2 ), (22) 

g=-oo 

where to, to' — 0, . . . , p — 1 and J m (/x, x c ) = (27r) _1 J Q 27r exp[im(a; c — u) + ipL cos(u)]du. It 
is clear from Eq. (22) that the matrix (21), and therefore the QE spectrum, is 27r-periodic 
in h at fixed /i. In addition, the matrix (21) is invariant under the "reflection" h — > 2ir — h 
at fixed \i provided x c and w\ satisfy the two relations 

x c TW\ = -x c ± wi + 27rj ± (23) 

for some integers j±. The solution of Eqs. (23) for x c is x c = + j-)/2. Since x c 
can be restricted to [0, ir/2] (see Sec. A), we see that the QE spectrum is invariant under 
h — > 27r — ft at fixed /i only if x c = 0, 7r/2. 

C. General quantum evolution in terms of periodic- wave-packet evolutions 

Following Ref. [4], let us express the momentum v in the form v = (n + f3)H, where n 
and P are, respectively, the integer and fractional parts of v/H. The (v, x c ) representation 
of a general 2D wave packet |$) for the system (1) will be then denoted by (p(n + (3, x c ). 
From the relation between the (u, x c ) and the (v, x c ) representations, one finds that 

$(w, x c ) — — J ip(v/H, x c ) exp(iuv/h)dv = J d(3exY>{i(3u)Qp{u, x c ), (24) 
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where $(«, x c ) is the (u, x c ) representation of |<3>) and 

00 

$p(u, x c ) — (f{n + (3, x c )exp(inu). (25) 

n=— 00 

Clearly, ^(w, x c ) is 27r-periodic in -u, so that Rel. (24) provides the decomposition of 
$(«, x c ) into Bloch functions exp(i(3u)$p(u, x c ) with quasimomentum f3. Now, since the 
evolution operator U (u, v; x c ) in Eq. (3) is 27r-periodic in both u and v, its application to 
such a Bloch function must "conserve" (3; in fact, it is easy to see that 

U(u, v; x c )[exp(iPu)^f3(u, x c )] = exp(i/3u)$^(u, x c ), (26) 

where 

$^(w, x c ) = U = u, v = h(3- ift^; a; c ^ ^(u, x c ). (27) 

The right-hand side of Eq. (26) is also a Bloch function with quasimomentum (3, since 
$^(w, x c ) is 27r-periodic in u due to Eq. (27). Rel. (24) then implies that the evolution 
of a general wave packet $(w, x c ) under U(x c ) can be decomposed or "fibrated" [4] into 
independent evolutions at fixed f3, 

U(x c )<S>(u, x c ) = I d/3esxp{i/3u)Up{x c )^p{u, x c ), (28) 
Jo 

where Up(x c ) is the operator appearing on the right-hand side of Eq. (27). Rel. (28) shows 
that the quantum dynamics (10) can be fully reproduced from that of periodic wave packets 
(25) under the corresponding operators Up(x c ). 

We denote by ip(v/h, x c ; t) = ^(x^ipi^v/h, x c ) the wave packet at time t in the (v, x c ) 
representation and restrict our attention, as in Sec. IV C, to the case in which one has a 
separation of the (v, x c ) variables, ip(v/H, x c ) = x ( x c)i>(v/fr), with normalized x an d i J - 
The expectation value {v 2 ) t of v 2 in ip(v/h, x c ; t) can then be expressed as follows: 

(v 2 ) t = j v 2 mv/h, x c ; t)\ 2 dv = J d(3Np (v 2 )^ t , (29) 
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where ip(v/h, x c ; t) = U t (x c )ip(v/h), 

fc2 00 

<*\* = AT E (n + /5) 2 |^(n + A x c ; t)| 2 , 



#„= £ |^(n + /3)| 2 , 



(30) 



and ip(n + /3, x c ; t) = U^(x c )^(n + (3). We notice that the expressions (29) and (30) are 
independent of x( x c) an d their x c -dependence is completely due to that of U(x c ) or Up{x c ). 
Fully analogous expressions can be written for the expectation values (u 2 ) t and (w 2 ) 7 v where 
7 is the quasimomentum characterizing Bloch functions 27r-quasiperiodic in the v direction. 
Let us now assume approximate quantum diffusions of (v 2 ) /3t and (u 2 ) yt for all (f3, x c ), 
(7, x c ), and for sufficiently large t, i.e., {v 2 ) /3t ~ 2Dp(x c )t and (u 2 ) lt ~ 2D ur/ (x c )t. The 
expectation value of the kinetic energy then exhibits the approximate diffusive behavior 
(lo 2 (u 2 + v 2 )/2) t « 2D K e(^ c )^ where 



Dke(xc) = (cu 2 /2) 



yo ^0 



Following the same reasoning as in the case of f3 — considered in Sec. IV C, the diffusion 

coefficients Dp(x c ), D U;1 (x c ), and -Dke(^c) for sufficiently small K are expected to satisfy 

approximately formula (18) with replaced by corresponding proportionality constants 
n( ± ) n(±) anr i 



(±) 

H,KE 



(^ 2 /2) 



1 d^Dg> + d-yN^Dg^ 
./o 



D. Cases in which D 



(+) 

H./9 
(+) 



We show here that there are cases in which = \D^^\ (similarly, one can show 



that there are cases in which = \). The constants D^L and \D^L\ are 

the diffusion coefficients for ('5 2 ) /3t under the standard Harper evolution operator Un — 
exp {i [cos(w) + cos(t))] /h} and its inverse U^ 1 , respectively. In the u representation, 



(-) 



>(+) 



»(-)i 



Uft(u) = exp 



cos(w) + cos ( ih 



d_ 

du 



(31) 



Consider the Bloch state A ( g(u) = exp(i/3u)$(u), where $(«) is given by the right-hand side 
of Eq. (15), and define the variable v! = 2u — u. It is easy to see that 
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Kp{u) = CA^u'), 



(32) 



where C = exp[2iu(v + f3)/h] is a constant phase factor. Choosing now u = jir, j integer, 
and using Rels. (31) and (32), we find that 



A (u; -t) = U^{u)K p {u) = C U^A^u') = CAUu'; t) 



(33) 



It follows from Rel. (33) that 



<« 2 > 



P-t 



-h 2 



2tt 



d 



(34) 



since (v 2 ) /3t must be real (and positive). The equality D^p = \D^p\ * s an immediate con- 



sequence of Rel. (34). 
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Table I. D(0)/K and D(n/2)/K 2 for several values of K. 



K 


0.15 


0.157 


0.17 


0.18 


0.19 


0.2 


D(0)/K 


0.458 


0.459 


0.461 


0.461 


0.460 


0.459 


D{ix/2)/K 2 


0.219 


0.226 


0.237 


0.243 


0.244 


0.243 
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FIGURES 



FIG. 1. Classical stochastic webs for K = 0.157 and a = tt/2, generated by the map (8) 
in three cases: (a) x c = 0; (b) x c = 1.47; (c) x c = ir/2. The range of both u and v in all the 
plots is [— 2n, 2n}. 

FIG. 2. "Hofstadter butterflies" (HBs) E(h) for fj, = 0.5 and: (a) x c = 0, -4 < E < 4; 
(b) x c = 1.47, -1<E<1; (c) x c = tt/2, -1 < E < 1. 

FIG. 3. Perturbed double HB #'(/*) (-2 < E' < 2), obtained by locally scaling the 
HB in Fig. 2(c) by e _1 = \pi sin(/i/2)] _1 . The perturbation is reflected by the fact that 
the spectral width AE'(ir) (at the HB center) is slightly smaller than the maximal width 
A£'(2tt) = 4. 

FIG. 4. Solid lines: (v 2 ) t for h/(2n) = [51 + (V5- 1)/2]-\ K = 0.157, and a; c = (upper 
line), x c = 1.47 (middle line), x c = ir/2 (lower line). Dashed lines: Corresponding classical 
quantity (vf) for the same values of x c (see text for more details). For the sake of visibility, 
the lines for x c = tt/2 were shifted below by multiplying both (v 2 ) t and (v 2 ) by 0.01975. 
The inset shows the normal plot of (v 2 ) t for x c = n/2. The values of H and K above were 
used also to obtain the results in Figs. 5 and 6. 

FIG. 5. Diamonds: Numerical results for the quantum-diffusion coefficient D(x c ) (see 
text for more details). Solid line: Least-square fit of formula (18) to the numerical data, 
£>Jj +) = 0.224. 

FIG. 6. Middle line: Global quantum diffusion of (v 2 ) t , given by a uniform average of 
(v 2 ) t over x c = kn/AO, k = 0, 1, . . . , 18 (see text for more details). The upper and lower 
lines correspond to {v 2 ) t for x c = (k = 0) and x c = 2n/5 (k = 16), respectively. 
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